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1.  Introduction 


The  exclusion  process  was  introduced  by  Spitzer  in  /lO/. 
Comprehensive  treatments  of  this  interacting  particle 
system  are  given  in  /7 /  and  /!/. 

In  the  exclusion  model  the  particles  attempt  to  move 
independently  according  to  a  Markov  kernel  on  a  given 
countable  set  of  sites.  But  any  jump  which  would  take  a 
particle  to  an  already  occupied  site  is  suppressed.  That 
means  there  is  always  at  most  one  particle  per  site. 

This  paper  is  a  continuation  of  /8/  where  we  considered 
the  so  called  wide  range  birth  and  death  exclusion  process 
in  random  medium.  We  proved  in  /8/  a  law  of  large  numbers 
for  this  measure  valued  process  and  derived  a 
deterministic  macroscopic  equation  describing  the 
evolution  of  the  occupation  rate  of  the  sites.  This 
macroscopic  equation  holds  under  rather  general 
conditions.  It  can  be  specified  by  the  choice  of  the 
Markov  kernel,  this  means  by  the  given  jump  intensity  of 
the  particles. 

Within  this  paper  we  will  study  the  case  of  a  nonsymmetric 
local  jump  intensity  which  allows  only  jumps  into  the 
neighborhood.  For  instance,  such  jump  intensity  is  of 
Interest  In  modelling  of  stochastic  charge  transport  (see 
/9/)  . 

In  the  following  we  will  show  under  suitable  assumptions 
that  the  asymptotics  with  respect  to  vanishing  mean  jump 
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size  yields  a  nonlinear  second  order  partial  differential 
equation  characterizing  the  evolution  of  the  limiting 
particle  concentration.  This  equation  represents  a 
continuity  equation,  where  physical  quantities  as  current 
density  vector,  conductivity,  static  potential  and 
chemical  potential  can  be  easily  identified.  It  seems  that 
this  continuity  equation  models  important  transport 
processes  in  physics,  chemistry,  biology,  electronics, 
social  sciences  and  other  fields. 

The  first  part  of  the  paper  generalizes  the  birth  and 
death  exclusion  process  introduced  in  /8/.  It  formulates  a 
law  of  large  numbers  for  the  case  that  the  mean  number  of 
sites  per  unit  volume  tends  to  infinity.  The  second  part 
derives  the  evolution  equation  for  the  limiting  particle 
concentration. 

2.  Exclusion  process 
2.1.  Random  medium 

In  the  following  we  will  introduce  a  generalization  of  the 
wide  range  birth  and  death  exclusion  process  considered  in 
/8/. 

For  unexplained  notations  and  definitions  we  refer  to  /4/ 
or  /5/. 

Let  iiQfFfP)  denote  the  basic  complete  probability  space. 
F  =  increasing  right-continuous  family  of 


complete  sub— <5  —■fields  of  F. 

Further,  B(E)  represents  the  Borel-cJ  -algebra  of  a 
topological  space  E. 

The  sites  are  located  in  a  closed  bounded  domain  Per  R  , 

de  ■{l,2,...  J  • 

We  introduce  a  finite  ^  —additive  measure  A.  on  B(Q)  which 
is  called  intensity  measure  of  sites. 

(N  )  denotes  a  sequence  of  -measurable  simple 

n  n  ^ 't 

counting  measures  on  B(Q) ,  the  so  called  counting  measures 
of  sites. 

The  parameter 

(2.1.1)  n=E  (0) <  j  dq) 

Q 

can  be  interpreted  as  the  mean  number  of  sites  per  unit 
volume. 

For  K  €  <0,  cx>  )  we  denote  by  the  set  of  bounded 

Lipschi tz-continuous  functions  f|Q  ->  [-K,K]with 

(2.1.2)  |f(u)-f(q)|  «=  K|q-u  1 

for  all  u,q  €  using  the  usual  Euclidean  norm.  Now  we 

assume  that  for  each  KC(0,oo)  it  holds 

(2.1.3)  lim  E  sup  (/  f  (q)  (^  -A  )  (dq)  ) 

n-><»  f  Q 

and  for  all  n>l  we  have 
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(2.1.4)  E(^  N  (Q)  )'^  =  K  <  oo 

n  n  O' 

The  counting  measure  represents  the  random  medium 

within  our  microscopic  stochastic  model.  The  above 

conditions  are  satisfied  for  a  wide  class  of  regular 
lattices  and  other  point  processes  including  the  Poisson 
point  process  (see  /8/). 

2.2.  Harkov ian  jump  machan ism 

We  denote  by  ^  the  counting  measure  of  particles  at 
time  t^O. 

Further,  we  introduce  an _F-adapted  cadlag  Poisson  counting 
measureyu^  on  B(CO,«®))®  B<S>  •  B(Q)  which  is  characterised 
by  its  dual  predictable  projection  which  is  here  its 
intensity  measure 

(2.2.1)  (dr  ,du,dq)  =n'^  w^  (u , q) (du) (dq)  dt . 

We  assume  that  the  jump 'rate  (u,q)  is  nonnegative  and 
Lipschi tz-continuous  with  respect  to  u  and  q  uniformely 
with  respect  to  t. 

The  counting  measure  generates  the  so  called  possible 

jumps  of  particles.  But  only  jumps  from  occupied  into 
vacant  sites  will  take  place.  Therefore,  we  will  have  at 
most  one  particle  at  each  site. 

2.3.  Birth  and  death  of  particles 

We  include  in  our  microscopic  stochastic  model  also  the 
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effects  of  birth  and  death  of  particles. 

The  possible  birth  ot  death,  resp. ,  of  particles  is 
generated  by  the  F-adapted  cadlag  Poisson  counting 
measures and^^  ,  resp.,  which  are  defined  on  B(  [O,**®))® 
g(Q)  and  characterized  by  their  dual  predictable 
projections 

(2.3.1)  T  <dt,du)  *  ;5\  (u)  N„<du)dt 
and 

(2.3.2)  V  (dt,du)  =  w^  (u)  N„(du)dt 

resp.  The  birth  rate  (u)  and  death  rate  (u)  are 

assumed  to  be  nonnegative  bounded  and  Li pschi tz -continuous 
with  respect  to  uf  Q  uniformely  with  resprect  to  t. 
Furthermore,  we  suppose  that andj^^  are  independent. 

A  birth  (death,  resp.,)  at  u  at  time  t  will  take  place 
only  if  u  is  vacant  (occupied,  resp.)  at  this  time.  So 
also  in  the  case  of  a  bir*th  or  death  it  is  ensured  that  we 
have  at  most  one  particle  at  each  site. 

2.4.  Initial  condition 

The  function  Y|Q  ->  [0,l]  denotes  the  initial  occupation 
rate  of  the  sites.  We  assume  that  at  time  t=0  at  most  one 
particle  is  at  each  site,  L is  F„  -measurable  and  for 

’  f),0  =0 

all  K  c  (0,oo  )  it  holds 

lim  E  sup  (/ f  (q)  (^  -  y’(q)A  )  (dq)  )^ 

n->oo  Q 


(2.4.1) 


“0. 


2,5,  Stochastic  equation 

Let  denote  the  Dirac  measure  at  u  i  Q  and  the  le-ft 

hand  limit  of  the  measure  of  particles  at  time  t>0. 

Now  we  can  define  the  measure  valued  exclusion  process 

=  {(-„  -i  unique  solution  of  the  following 

stochastic  equation  (see  /6/) 


<  {ujy  il-L^  ^  (  {q}  )  (ds,du,dq) 

t 


^  ^  y^„(ds,dq) 

o  Q  y  ' 

<fu}yjd6,du), 

which  describes  the  evolution  of  L_  driven  by  jut  , 
and  .  One  easily  notes  how  the  ^’-adapted  cadlag 
piecewise  constant  and  Markovian  measure  valued  process  L, 
remains  at  any  time  t  with  at  most  one  particle  at  each 
site.  Furthermore,  the  interaction  between  the  particles 
caused  by  the  exclusion  mechanism  is  reflected  by  the 
logistic  nonlinearity  L  <1-L  )  within  the  second  term  of 
the  right  hand  side  of  equation  (2.5.1). 

2.6,  Occupation  rate 

Let  Q  denote  the  support  of  the  intensity  measure  of 
si tes/\..  Now ,  for  t-0  and  qe  Q  we  introduce  the  so  called 


occupation  rate  H(t,q)  which  we  will  later  interpret  as 
the  probability  that  a  site  at  q<Q  is  occupied  at  time  t. 
We  define  the  occupation  rate  H  as  the  unique  solution 


<see  /8/)  of  the  following  integro-di f f erenti al  equation 

<2.6.1)  ^  H<t,q)=»<l-H(t,q)  </w^(u,q)H(t,u)  A  (du)+w^(q)  ) 

f  S 

-  H(t ,q) (y  <q,u) (1-H(t,u) ) 4  (du)+w^ (q) > 

3 

for  all  t>0,  qfQ,  with  initial  condition: 

(2.6.2)  H<0,q)=  'j  (q)  , 

for  all  qfQ.  It  can  be  shewn  as  in  /6/  that  we  have  for  all 
t^O  and  qfQ 

(2.6.3)  H<t ,q)  C  [O, 1 ]  . 

Equation  (2.6.1)  is  a  bi lance  equation  for  the  macroscopic 
evolution  of  the  occupation  rate  and  allows  the  following 
interpretation:  The  occupation  rate  changes  its  value  in 

dependence  on  the  occupation  rate  at  other  points.  It 
increases  at  qf  Q  proportional  to  the  non-occupation  rate 
(1-H(t,q)>  and  the  sum  pf  the  birth  rate  ^^<q)  together 
with  the  occupation  rates  H(t,u)  for  the  sites  at  ufO\  fq} 
which  are  weighted  by  the  intensity  w^  (u,q)  for  jumps  from 
these  sites  into  q.  On  the  other  hand  H(t,q)  decreases 
proportional  to  its  own  actual  value  and  the  sum  of  the 
death  rate  <q)  together  with  the  non-occupation  rates 
(1-H(t,u)  of  the  other  sites  at  uf  5\  {q}  which  are 
weighted  by  the  intensity  w^ (biU)  for  the  jumps  from  q 


into  these  sites. 
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We  remark,  that  Qroeger  proved  in  /3/  that  a  unique 

solution  exist  for  the  bilance  equation  (2.6.1)  together 
with  the  Poisson  equation  describing  the  self  consistent 
static  potential . 

2.7.  L»m  of  large  numbers 

Now,  we  can  formulate  a  law  of  large  numbers  for  the  above 
introduced  wide  range  exclusion  process  in  random  medium 
for  the  case  that  the  mean  number  n  of  sites  per  unit 
volume  tends  to  infinity,  n  ->  oo  . 

Under  the  above  assumptions  it  holds  for  fixed  T,K  e  (0,<») 
(2.7.2)  lim  E  sup  E^sup  (/f  (q)  (-L  -H  ( t ,  q)/t )  (dq)  )*  I  fJ=0 

n->oc,  ft  O^t^T  Q 

The  choice  of  the  class  of  test  functions  C  and  the 

K 

positioning  of  the  expectations  are  cruical  for  the  proof 
of  the  theorem  which  uses  semimartingale  methods.  It  can 
be  omitted  because  it  would  be  almost  the  same  as  that 
which  is  given  in  /&/ . 

The  above  law  of  large  numbers  shows  that  the  random 
functional 

'-..I 

Si 

converges  in  the  mean  square  sense  uniformely  with  respect 
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to  the  time  te[0,Tj  and  all  test  functions  ‘f  (  C ^  for  n  ->  «? 
to  the  deterministic  functional 

J'f  (q)H(t  ,q)/t  (dq)  . 

The  inner  conditional  expectation  in  (2.7.2)  relates  to 
the  driving  Poisson  counting  measures  and^^^.  The 

remaining  outer  expectation  averages  the  random  medium  and 
the  random  initial  occupation  of  the  sites. 

3.  Asymptot ics  of  the  mmcroscopic  nonequx I ibr ium  dynamics 

3.1.  Specifications  for  the  case  of  a  local  jump  intensity 

For  simplicity  let  us  cho  se  the  domain 

(3.1.1)  Q=(0,  ])>:..  .X  (0,1-  ) 

1  « 

wi  th 

(3.  1.2)  1*  ^  (0,oo  ) 

for  all  i  f  {l...,dj  . 

QQ  is  the  boundary  of  Q  and  Q  denotes  its  closure. 
Further,  ? Q'  is  that  part  of  3 Q  which  does  not  contain 
corners  or  edges. 

For  i,jf  we  denote  by  the  set  of  all  i-times 

continuous  differentiable  functions  on  Q,  and  3  is  the 
set  of  all  functions  on  ^  XD  which  are  i-times 

continuous  differentiable  with  respect  to  the  first  and  j- 
times  with  respect  to  the  other  coordinates,  T  (  (O,  ). 


-  1 1  - 

We  denote  by  the  outward  unit  normal  of  t?  q'  at  q  f7)Q,  by 
div  a  the  divergence  of  a  vector  a,  by  a-e  the  usual 
scalar  product  between  two  vectors  a  and  e  and  by  grad  f 
the  gradient  of  a  function  f  on  £ 

We  assume  that  the  intensity  measure  of  sites  is 

absolutely  continuous  and  we  write  for  all  qf  Q. 

(3.1.3)  (dq)  =X(q)dq. 

We  suppose  that  there  exists  a  version  of  the 
concentration  of  sites^  with  the  properties 


(3.1.4)  Xiq)  ^  >0 

for  all  q  <■  Q,  where 

(3.1.5)  ^  e  c'’  . 


Under  this  assumption  we  note  that  the  support  Q  of  A 
coincides  with  To  formulate  the  local  jump  intensity  we 
introduce  the  Lipschi tz- continuous  probability  density 


(3.1.6) 


p(q)=exp  <■-  q  } 


(  f  exp  {  -  |u|j  du) 


for  all  q  f  . 

We  remark  that  the  proposed  approach  would  also  work  for 
other  probability  densities  if  they  show  moment  properties 
as  those  listed  in  Section  4.1. 


For  each  value  of  a  parameter  (0, 1 )  we  specify  the  local 
jump  intensity  for  all  t-0  and  u,qf  Q  by  the  expression 
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(3.1.7)  w  (u,q)  =  Tr<t,q)^  TrCt.u)^  p  <-  (u-q)  )  << 

t  ^ 

wi  th 

(3.1.8)  'rr(t,q)  =  e>:pj  -^(t,q)J  , 

where  ^  (t,q)  is  the  so  called  static  potential  and  we 
assume 

(3.1.9)  P  €  C"'"'  . 

We  note  that  -for  smaller  the  Jump  intensity  w^  is  more 
local i zed. 

To  interpret  the  jump  intensity  w^  we  remark  that  one  can 
show  by  the  use  of  the  properties  of  the  probability 
density  p  listed  in  Section  4.1  that  we  have  for  q«  Q  the 
asymptotic  drift  vector 

lim  y(u-q)w  (q,u)du  =  -b  grad^(t,q) 
cC  ->0  0 

and  the  asymptotic  variation 

lim  J  (u-q)  w. (q,u)du  =bd, 

aC->OQ 

where  b  is  a  positive  constant  depending  on  d.  For 
instance  we  have  for  d=l  the  value  b=l/2  and  for  d=3  the 
value  b=T  . 

Now,  the  occupation  rate  depends  on  the  parameter  oC  c  (0, 1 ) 
and  we  use  in  the  following  also  the  notation 
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(3.1.6)  H, <t ,q)=H(t,q) 

0C 

for  all  tf  [  0,t3  and  q  €  Q • 

3.2.  An  integral  quation 

In  the  following  we  characterize  the  limiting  dynamics  of 
the  occupation  rate  for  For  this  purpose  we 

assume  that  converges  pointwise  to  a  function 

H  |[0,t3XQ->  [0,1]  such  that 

(3.2.1)  lim  (H  ,  (t,q>-H(t,q>>=0 

->0 

for  all  tf  [OjT]  and  q  f  Q. 

Further,  we  assume  that  the  time  derivative  of  the 
occupation  rate  is  uniformly  bounded  for  alleCf(0,l), 
t  f  C  0,T]  and  q  f  Q  wi 

(3.2.2)  (t,q)  I  . 

Finally,  let  us  suppose  that  there  eKists  a  constant 

(3.2.3)  »f€(0,l/2) 

such  that  for  al  1  «(.  C  (0, 1 )  ,  qfQ  and  f  [0,T]  we  have 

(3.2.4)  H^(t,q)f  [af,l-d€]. 

That  means  the  occupation  rate  is  never  O  or  1. 


^  i  *1 

Let  us  denote  by  ^  the  set  a-f  -functions  f  f  C  with 
<3.2.5)  ^  ♦<q.=0 
fo  all  q  e  VQ*. 

Now,  we  are  able  to  characterize  H  as  solution  of  an 
interesting  integral  equation. 


The  limit  H  of  the  occupation  rate  for«< ->0  satisfies 
under  the  above  assumptions  for  all  fe  and  tf  C 0,Tj  the 
integral  equation 


(3.2.7)  /f  <q)  <H<t,q)-H<0,q»X  (q)dq 
Q  ^ 

<=  jf{\  H  <B,  q)  [  di  V  ( .^  (q)grad  f  (q)} 

-^^(q)  ( 1-H(s,q)  )grad  f  (qlgrad  ^  <s,q)J 

+f  (q)  [  <l-H(B,qJ^  (q)-H(B,q)  q>J^  <tl)J  db  ds- 


The  proof  of  this  theorem  is  given  in  Section  4.3. 

Equation  (3.2.7)  gives  a  characterization  of  the  limit  H 
which  avoids  smoothness  assumptions  on  H.  Therefore  one 
can  say  that  (3.2.7)  gives  a  rather  weak  description  of  H. 
Under  sufficient  smoothness  assumptions  we  will  show 
within  the  ne>st  section  that  H  is  the  solution  of  a 
corresponding  nonlinear  partial  differential  equation 


which  allows  a  direct  interpretation  of  the  dynamics 
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already  described  by  the  equation  (3.2.7). 

3.3.  Asymptot ic  occupation  rate 
We  assume  the  initial  condition 

(3.3. 1)  H(0, . )  f  C’ 
wi  th 

(3.3.2)  R(o,q)+TT(o,q)  (l-Tr(o,q>  );^  ^  (o,q)=0. 

J  ■?  ^ 

%  ? 

Now,  we  are  going  to  introduce  a  -function  R  on  xQ 

which  we  call  asymptotic  occupation  rate.  We  suppose  that 
there  exists  a  function  R  f  £.  which  is  the  unique 
solution  of  the  nonlinear  partial  differential  equation 


(3. 

3.3) 

R<t, 

q)=:  X^(q)^div<  .^Nq) 

(grad 

R(t ,q) 

+R(t , 

q)  ( 1-R (t  ,  q)  > 

grad  j 

^ (t ,q) ) > 

-Kl-R 

(t  ,q)  I'w  ^  (q) 

-R(t , 

q)w  ^  (q) 

for 

all 

t  >0  and  q  £  Q 

with  reflect 

ing  boundary  condition 

(3. 

3.4) 

0— 

,q)+R<t 

,q) (l-R<t ,q) 

f (t ,q)=0 

for 

all 

t-0  and  q  £  Oq 

* ,  and  initi 

al  condition 

(3. 

3.5) 

R  (0,  q) 

=H(0,q) 

for 

all 

qf  Q. 

We  note  that  (3.3.3)  represents  a  generalization  of 
Burger's  equation. 
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The  following  theorem  shows  under  sufficient  smoothness 
assumptions  on  the  limit  H  of  the  occupation  rate  for 

->  0,  that  Tl  coincides  with  the  asymptotic  occupation 

rate  R. 

ItieoC§!D 

If  we  assume  the  property 

<3.3.7)  Hf 

then  we  have  for  all  <t,q)f  ]  >’>0  ihe  e  v' «  ie  n  c€“ 

<3.3.8)  H<t ,q)=R<t ,q) . 

The  proof  of  this  assertion  is  given  in  Section  4.4. 

We  remark  that  the  smoothness  assumption  <3.3.7)  on  TT 
could  be  considerably  weakened  by  an  appropriate 
functional  analytic  formulation  of  the  nonlinear  partial 
differential  equation  <3. 3. 3) - <3. 3. 5) .  For  instance  such  a 
weaker  formulation  could  be  based  on  methods  described  in 
/2/  or  / 6 /-  Here  we  have  chosen  rather  strong  smoothness 
assumptions  on  R  and  H  to  derive  the  dynamics  described  by 
the  equations  <3.3.3)  -  <3.3.5)  without  technical 

difficulties  in  the  formulation  and  prove  of  Theorem 
3.3.6.  It  remains  an  interesting  problem  to  derive  these 
equations  under  much  weaker  assumptions. 

Finally,  we  remark  that  the  result  could  be  generalized 
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to  the  case  of  a  nonisotropic  probability  density  p  and  a 
general  regular  domain  Q. 


3.4,  Continuity  equation 

For  better  interpretation  of  the  asymptotic  occupation 
rate  R  together  with  other  physical  quantities  as  particle 
concentration  and  current  density  vector  we  rewrite  the 
equation  (3.3.3)  in  the  form  of  a  continuity  equation. 

Let  us  introduce  for  all  t^O  and  q  f  Q  the  particle 
concentration 


(3.4.1)  ^(t  ,q)»R(t,q)  A.  (q) 


and  the  vector  function 


(3.4.2)  j  (t  ,q)=-^;^^(q)  (grad  R(t,q) 

+R(t  ,q)  (1-R(t  ,q)  )  grad^(t,q)) 


which  we  will  call  current  density  vector. 

Then  it  follows  from  (3.3.3)  the  continuity  equation 

-n 

(3.  4. 3)^  ^  (t ,  q)  =  -div(j(t,q)) 

+  (-^  (q)  -  ^  <t,q>  )'w'^(q)  -^(t »  q)  w  ^  (q) 
for  all  t>0  and  q f Q,  with  reflecting  boundary  condition 


(3.4.4)  J(t,q)*'?  =0 

% 

for  all  t-0  and  q^pQ',  and  initial  condition 
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(3.4.5)  ^  (0,q)=H(0,q)  ^ (q) 

for  all  qeQ. 

The  above  continuity  equation  relates  the  time  derivative 
of  the  particle  concentration  with  the  current  density 
vector . 

We  can  also  write  the  current  density  vector  in  the  form 

<3. 4. A)  i<t,q)» 

-  j  A<q)R(t,q)  <  1-R  (t  ,  q>  >  grad  <J(t  ,q)+  J  (t,q>  )  , 
for  all  t-0,  qf  Q,  where 

(3.4.7)  J(t,q)=ln(R(t,q>  (1-R(t,q>  > 

is  the  so  called  chemical  potential  at  t-0  and  qc  Q. 

One  notes  that  the  current  density  vector  shows  into  the 
opposite  direction  of  the  gradient  of  the  sum  of  the 
chemical  and  static  potential. 

The  length  of  the  current  density  vector  is  proportional 
to  the  so  called  conductivity 

(3.4.8)  <3(t,q)=  J  b^^(q)R(t,q)  (l-R(t,q)  ) 

at  t-0  and  q€  Q,  which  contains  a  logistic  nonlinearity 
with  respect  to  the  asymptotic  occupation  rate  R.  The 
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conduct! vi ty 

reaches 

its  max i mum  at 

an 

asymptot i c 

occupation  rate  with  value  R=l/2  which  corresponds  to  the 

value  of  the 

chemical 

potential  ^  =0. 

The 

logistic 

nonl i neari ty 

of  the 

conductivity  is 

caused 

by  the 

i nter action 

of  the 

particles  within 

the 

exc 1 usi on 

dynamics,  where  jumps  to  occupied  sites  are  excluded.  As 
important  result  o-f  this  nonlinearity  it  follows  that  the 
particle  concentration  is  bounded  by  the  concentration  of 
sites  and  we  have  for  all  t-0  and  Q 

(3.4.9)  0^  57<t  ,q)^  ^  (q)  . 

Finally,  we  note  that  there  is  no  current  through  the 
boundary  VQ'  what  is  natural  in  our  model.  It  seems  that 
the  continuity  equation  <3.4. 3) - (3. 4. 5>  can  be  used  to 
model  important  transport  processes  with  birth  and  death 
effects  in  physics,  chemistry,  biology,  social  sciences 
and  other  fields. 

4.  Proofs 

4.1  Properties  of  the  probability  density 

For  easier  reference  we  list  some  properties  of  the 
probabi 1 i ty  densi ty 

p(q)=exp  {  -  (  q  (]  (  exp  {  “  |  u  (  }  du> 

introduced  in  (3.1.4). 

cf 

p  is  symmetric  and  we  have  for  all  q€  R  that 
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(4.1.1)  p(q)=p(-q)) 


Far  all  i  ,  Jc  f  ^  1  , .  .  .  ,d  }  one  obtains  by  the  use  of  the 
abbreviation 

p=p  <u-q)  )  oT 

the  following  moment  properties 


(4.1.2)  lim  ^(u  -q-)p  du=0  for  qt  Q 

oC->0  Q  ^  ^ 

(4.1.3)  lim  /  (u-q-lp  du=d-  <q)  for  q^jQ, 

o<.->o  Q  ^ 


(4.1.4)  lim  J  (u ■ -qj  )  (Uj^ -q^^  ) p  du=  j 
oc~>o  0  ^  ' 

and 


if  q  (  Q  and  i =k 
if  q  f  Q  and  i /k 
(q)  if  q(  0  Q 


(4.1.5)  lim  /  |u-q|  p  du*  0  for  all  qt  O  and  lf{3,4,...} 

oc  -*0  5 

where  b  is  a  fixed  bounded  positive  constant  depending  on 


the  dimension  d  and  it  holds  for  d:  and  Cj ,  the  estimate 

l.K 

(4.1.6)  |c_j^  (q)|  +|dj  (q)  I 

for  al  1  q€t?  Q. 


4,2.  Some  preparations 

For  all  aC  €  (0,1),  tf  [  O,  T  ]]  and  qf  Q  we  define  the 
function 

(4.2.1)  A^(t  ,q)«=H^  (t  ,q)  (1-H^  (t,q)  )'’'7r  (t  ,q)'’  ffo,®©] 

and  it  follows  with  (3.2.4)  that 

(4.2.2)  H^(t,q)=A^  (t,q)-JT  (t,q)  (1+/^  (t,q)Tr  (t,q)  )“’  f  [^e,l-ae  ] 
Now,  we  get  from  (4.2.1)  with  the  estimate  in  (4.2.2)  and 
the  boundedness  of 'rr  because  of  (3.1.B)  and  (3.1.9)  also 


1 
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an  estimate  for  with 

14.2.3)  (t,q)  f  .ri^] 

and 


<4.2.4)  0<  '' 


for  all  o<-f(0,l),  tfrO,TJ  and  q  6  Q. 


With  the  notation  (4.2.1)  and  the  explicit  form  of  the 
local  jump  intensity  w^  in  (3.1.7)  we  can  write  the 
integro-di f f erenti al  equation  (2.6.1)  in  the  form 

(4.2.5)  ^  <t,q)=  /(A^  (t,u)-A^  (t ,  q)  )  B  A  (u)  du-H^ 

O 

for  all  o(f(0,l),  tf  <0,Tj  and  q  f  IT  with  the  abbreviations 


(4.2.6)  B=<l-H^(t,q))  (1-H^<t,u))  Tr(t,q)^  ■Fr<t,u)^  p 
and 

(4.2.7)  =(1-H^  (t,q)  )w'^  ( q) ( t ,  q)js  ^  <  q>  • 

To  simplify  our  natation  let  us  use  also  the  following 
abbrevi ations 

f^  =f  (q)  ,  =  A  (q)  ;  A^=A^(t,q),  =H^  (t ,  q)  ,  ■n^=H  (t ,  q)  , 

TTy  =7r(t,q)  ,  J  (t,q)  . 

Further,  we  suppress  the  indication  that  we  integrate  over 
the  set  Q. 

Now,  it  follows  from  (4.2.5)  and  the  symmetry  of  B  for  all 


L 


« 


6(0,1),  tt  <0,TJ  and  i  t  C  that 


(4.2.6)  Jf^  -A^  >B;i_du.-|^dj*/f,  F,  dt, 

'//  ”,  '£<«  dy. 

It  holds  the  following  continuity  property  for  A.  . 


Lemma  4^2^? 

There  exists  a  constant  eiO^oo)  such  that  for  each 
t  €  (OfTj  and  o6f  (0,1)  we  have 


(4.2.10)  (A..-A.)'  p  dudq  . 


^  u  ■  ■  J 


Proof 


If  we  set  in  (4.2.8)  f^  =A^  ,  then  we  obtain  with  (4.2.3), 

(3.2.2),  (4.2.7),  (3.2.4)  and  the  boundedness  of  vT  and  w, 

t  -i 

that 

(4.2.  11)  ff  (A^  -A  )  ^  B2  dudq=  -  f  A  +F  );?  do^K* . 

Because  of  (4.2.6),  (3.1.4),  (3.2.4),  (3.1.8)  and  (3.1.9) 

we  can  conclude 

(4.2.12)  ff  (A^-A^)^  P  dudq-K'(K^  K '  '  ^  .  O 

Using  this  result  we  obtain  also  an  asymptotic  continuity 
property  for  H  : 


Lemma  4.2. 13 


For  all  t  6  (0,TJ  we  have 


<4.2.14)  lim  Jl  ^  !  u  "Q  I  P  d'-id  q  ~0 . 

oC  -  >  o'-'  “ 

Proof 

l-y  the  U5>e  ot  (3.2.4),  (4.2.2),  (3.1.8),  I  en.ma  (4.2.9)  and 

(4.1.5)  it  follows  for  all  t  f  (0,Tjth3t 

(4.2.15)  lim^(H_^-H  )  |u  — q|  p  dudq 

o£  ->0  "*  ^ 


1  i  m  ^2  1  H  -H  1  ju-qj  p  dudq 
-;0 

lim^''2(A^i37^  -TT  1  +  IrjA^^-Aj) 

1 

|u-q  1  p  dudq 

=  lim  {/Tk'  I  q-u  pdudq+yThC"  /  A  -A^  |  I  u-q|^  pdudq 
^  ->0 

=  0  +  1  im  K  <  /  <  (u-q)^  p  du)  dqfdq)^ 

J.->0  '  '  t  J 

~  1  i  m  K  (  j{  Jj  ~^u/  P  ^  ( (u-q)^  p  du)  dq)^ 

-  >0  f  ' 

=  o .  D 

4.3.  Proof  of  Theorem  3.2.6 

We  obtain  from  (4.2.8)  with  (4.  2.  1)  for  all  o^f(0,l), 

t((o,T],  and  ff^  by  symmetry  arguments 

(4.3.1)  /f^ 

=//  -Vr 

=  /[<^  -^u  ’rf  <H,-H^)+(l-H  )H,  (Tt/tt:^  -1) 

“  <  1  (T  *■  P  r7  dudq 


=  I  ( -f  -f  )  ■  H  ■+  '  1  h  >  H  .  *  - 


)  -1  ) 


+  <  1  -M  )  H  ,  (  1  -sy.p  •  r  <  j  -  0 x  P  '7  A.  dudq 

-  y  ^  ^  -  ‘'I 


I  a  -f  )-/  -H  +  -S^)  L  (1-H 

•♦■  (  I  -  H  )  H  .  I  <■  R  '  p  ^  ^  d  LI  d  q 

J  ■'if. 


'i  W  H  +  -if  )[  H  (1-H  ) 

+  (H_^-H^)J  +R  J  dudq 


(1-H  )+H  (1-H^) 

K  ^  ^ 


(f.  -f,  ) 


{ 

R  )  )  p  dudq 


wi  th 


(4.3.2)  ^K^/u-q). 

How,  we  expand  by  the  use  of  the  Taylor  formula  the  given 

functions  f,  ^  and  ^  ,  apply  Lemma  4.2.13  and  use  the 

property  (4.1.5)  that  higher  moments  than  the  second 

(Tioment  of  p  vanish  asymptotically.  Then  it  follows  from 

the  abo'/e  result  with  assumption  (3.2.1)  and  the  other 

~  2 

fTioment  properties  of  p  for  all  t  (  [  O ,  ~f  ]  and  ^  ^  that 

(4.3.3)  I=|f  (Fr<t ,  q) -H(0,q)  )^  dq- //f  r  (  1  -H(s,  q)  )  w} /q) 

-H(s,q)w,-  (q)J.l  dqds 

f  t 

=  1  i  m  (  /  <H^(t,q)  -H^(0,q)  j  /  f^  ^^dqds) 


1  1  (T,  ,  .  ;  -  _  —  t  ( M  q  .1 

^  -  O'  1=1'  • 


1  , 


- -  -f  (u  -q,-  )  (u  -q , )  ) 

i -v  -  ^  ■  i  ^  ■*  J 

•>1  *] 


^  d  p  “ 

1  T<II 

^  i=l  c^fj  ^  ^ 


d  p 

<-^  ^  (u--q  ))p  dudqds. 

?  r^l  ^  ^ 

We  note  that  for  f  f  C  we  have  for  all  q  t  '  the  boundary 
condi t i on 


f  =0. 


PV^  ? 

Now,  by  the  use  of  the  moment  properties  of  p  ;  (4.1.2)- 

(4.1.5)  and  assumption  (3.2.1)  it  follows  from  (4.3.3)  for 

all  tfC0,TJ  and  f  f  ^  ^  the  assertion 
i 

(4-3.4)  1=  //b[  -(grad  f^)  j(gradJ^ 


-Kgrad  ^ (grad>?^  ) 

+  7(div  grad  f  )  dqds 

//rH^rdiv(^y  gr^*sd  f^  ) 


-  y  (1  -H^  )grad  f^  •  grad  J  dqds. 


O 


4.4.  Proof  of  Theorem  J.J,6 

, 

From  Theorem  3.2.6  we  obtain  for  al I  t  f  (0,T]  and  f  €  ^ 


the  equation 


•A 


(4.4.  W  V  (-!-_-  _■  M  'h,  )  w.  ‘.c-  W.  (q>J  ) /„cq 

=  ,  f  ^  d  1  V  (  g  r  a  d  f ,  ) 

-  1  -Rj  )  (gr ad  f._  )  •  (gr  ad/^  )  J  dq 

-  f  ...  f^'~  4-^. ) 

1-1 


r 


'Vi  ^ 


c;  0  *■  ‘  J 

-H  )  )  (^1^)]  dq^  ...dq^. 


and  it  follows  by  partial  integration 


y  d  p  1  c?  - 
-/■••/ x>—  V 


?i 


-•H  (l-H  )-:^$  )dq  .  .  .dq  , 

d  1*/  ii'.’  1  •.,  •!•»  ,  T  ,  ,  _  ■>/'.. 

i  =  l  0  CO  j  •>  ^  -' 


where  we  used  the  abbr evi at i ons 


and 


V-  —  (q  y...|qt  ^1’  C/i  yas.^q;] 

1  1  T7  '■  i  1  f  f  « 


V. -<q^, - ,q^  ^  ,0,q^-^^  ,...q^  ). 


Because  of  (3.2.5)  we  can  write  for  all  t  6  (0,Tj  and  f  f  C 


<4 


.4.2)  V=  J  ^  (qr^d  (grad  ( 1  -H^  )  gr  ad  ^^  )  dq. 


i 


? 


Again  by  partial  integration  we  get  from  (4.4.2)  together 
th  (4.4.1)  for  all  t  €  (0,tD  ^nd  ffC^  the  resul  t 


wi 


( ^  .  3 )  '  -f  (  K,  -  ( 1  -  H,  >  <  q )  -R^  w_  <  q  >  j )  R  d q 

\  '  '  *  _  _  '  *  ’ 
f  r  di  V  <  <  \gr  ad  +H^  ( 1  -H^  )  gr  ad  )  )  dq 

1  ’ 


c 


+  J  i  H  +H  (1-H  )-f7-  f  )dq. 

:  c‘  ^  ^  ^  ^  J  -V 


"7^7 

Because  the  above  equation  holds  for  all  ft  C  and  we  assu.T.ea 
that  H  is  from  it  follows  for  H  with  (3.1.4)  and 

(3.3.2)  the  partial  differential  equation 


(4.4.4) 


) 


=  ~  di V (  ^  (grad  H  +H  (1-H  ) grad  i  ) ) 
1  ?  7  1  %  “% 


for  all  qt  Q  and  t  f  (0,TJ  with  the  boundary  condition 

for  all  qfQO'  and  t  f  f  0,T  3  • 

Now,  we  can  conclude  from  (4.4.4)  and  (4.4.5)  that  H  is 
the  unique  solution  of  the  equations  (3. 3. 3) - (3. 3. 5)  which 
means  for  all  t  (  (O,  T  3  and  q  ^  Q  that 


(4.4.6) 


H(t,q)^R(t,q)  . 
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